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Abstract 

In the nonlinear 0(N) sigma model at N = 3 unexpected cutoff 
effects have been found before with standard discretizations and lattice 
spacings. Here the situation is analyzed further employing additional 
data for the step scaling function of the finite volume massgap at 
N = 3, 4, 8 and a large iV-study of the leading as well as next-to- 
leading terms in 1/N. The latter exact results are demonstrated to 
follow Symanzik's form of the asymptotic cutoff dependence. At the 
same time, when fuzzed with artificial statistical errors and then fitted 
like the Monte Carlo results, a picture similar to N = 3 emerges. 
We hence cannot conclude a truly anomalous cutoff dependence but 
only relatively large cutoff effects, where the logarithmic component is 
important. Their size shrinks at larger N, but the structure remains 
similar. The large N results are particularly interesting as we here 
have exact nonperturbative control over an asymptotically free model 
both in the continuum limit and on the lattice. 



1 Introduction 



The 0(N) nonlinear sigma models in two space time dimensions are a very 
precisely studied class of quantum field theories in many respects. For N > 3 
they are asymptotically free like QCD and hence furnish a simplified test 
laboratory for all kinds of questions. This includes the question of the de- 
pendence of physical results on the lattice spacing and this is the subject 
studied in this publication. 

While we are interested in results at zero lattice spacing, nonperturbative 
Monte Carlo results are always at finite a. One hence has to assume an ana- 
lytical form for the asymptotic a-dependence to extrapolate to zero. Ideally 
this behavior is checked with good precision for a large range of a and then 
assumed for the 'remaining' extrapolation. In principle a systematic error 
has to be estimated here, which however is often deemed to be negligible 
compared with statistical errors, which tend to be amplified in the extrap- 
olation. This balance depends on the 'extrapolation fit' one uses. If it is 
restrictive with few parameters (but compatible with the data) statistical er- 
rors remain small but the danger of systematic effects is high. In the opposite 
case one may 'give away' statistical precision and overestimate total errors. 
It is presumably clear from these words that in the extrapolation step there is 
some danger of subjective elements getting into a presumed 'first principles' 
calculation. The situation is aggravated in QCD with dynamical fermions, 
because due to algorithmic and CPU-power limitations one cannot vary a 
over a large enough range. Continuum results then mainly arise under the 
assumption of a certain a-dependence. 

Almost universally in lattice field theory the model of the asymptotic 
cutoff dependence is derived from Symanzik's analysis of lattice perturba- 
tion theory to all orders [1J. It extends the investigation of renormalizability 
to subleading orders in the cutoff and is assumed to hold in structure also 
beyond perturbation theory and to thus ascertain the existence of the con- 
tinuum limit and predict the rate at which it is approached. The results for 
the rate are powers in a (usually a 2 or linear for fermions) modulated by only 
slowly varying logarithmic factors. 

In view of this situation there have been efforts to explore this in detail 
in the sigma model lab. The situation there is vastly better than in QCD: 
No critical slowing down due to cluster algorithms [2J , well-understood and 
well-measurable finite size quantities [3], reduced variance estimators [4 . The 
Symanzik prediction for this bosonic theory is clearly a 2 scaling. In earlier 
investigations [HI Ej it has been found that some results of the 0(3) model 
do not fit very well with this standard expectation down to rather small a- 
values. These studies are extended here in the direction of larger iV-values, 
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Figure 1: Lattice step scaling function E(w , N,a/L) normalized by the continuum value 
a(u ,N) at M = 1.0595. 



since it is known analytically that eventually at N = oo the Symanzik picture 
holds. We here report in Sect. 2 Monte Carlo results for N = 3,4,8 with 
the standard lattice discretization. In Sect. 3 we investigate the large N 
expansion including the leading correction of order 1/N. In Sect. 4 we draw 
some conclusions. 



2 Monte Carlo Data 

We consider the lattice O(N) nonlinear a- model 

Z= f l[d N s5(s 2 -l)e- s ^ (1) 

with the standard lattice action 

5 ( s ) = ^E(<v) 2 ( 2 ) 

where is the forward difference operator and 7 the rescaled coupling (fixed 
at large N). Both in our MC simulations and the large N expansion we 
determine the step scaling function (for a factor two rescaling) 

E(u,N,a/L) = M(2L)2L, u = M(L)L, (3) 
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where M(L) is the mass-gap [7j of the transfer matrix. It is determined from 
the asymptotic exponential fall-off of the two-point function at zero spatial 
momentum. In principle this refers to a periodic box with finite spatial 
size L and infinite extent in the time direction. In the simulations this was 
approximated as described in [3]. The renormalized coupling M(L)L used 
here differs by a factor (N — l)/2 from g 2 {L) in [3] to keep it finite in the 
large N limit. The lattice step scaling function is expected to reach a finite 
continuum limit, 

lim E(m, N, a/L) = a(u, N) (4) 

and we discuss here how this is approached. At large iV we shall use the 
expansion 

E(it, N, a/L) = E (u, a/L) + -J- E^u, a/L) + 0{N- 2 ). (5) 



L/a 


E(w , 3, 


E(u ,4,a/L) 


E(ii , 8, a/L) 


Y, (u ,a/L) 


Ei(n ,a/L) 


5 


1.29379(8) 
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1.31256(47) 


1.35065(40) 


1.379491 


-0.194123 


8 




1.31068(43) 


1.34891(35) 


1.378215 


-0.196201 


10 


1.27994(9) 


1.30825(39) 


1.34835(34) 


1.377524 


-0.197668 


12 


1.27668(9) 


1.30608(37) 


1.34721(32) 


1.377112 


-0.198665 


16 


1.27228(12) 


1.30477(34) 


1.34725(30) 


1.376665 


-0.199867 


24 


1.26817(9) 


1.30379(32) 


1.34649(28) 


1.376309 


-0.200939 


32 


1.26591(9) 


1.30263(31) 


1.34635(27) 


1.376170 


-0.201394 


64 


1.26306(16) 


1.30225(32) 


1.34570(25) 


1.376022 


-0.201921 


96 


1.26216(18) 


1.30138(33) 


1.34577(24) 


1.375990 


-0.202042 


128 


1.26187(13) 


1.30149(34) 


1.34563(24) 


1.375978 


-0.202089 


cont. 


1.261208(1) 


1.3012876(1) 


1.345757(2) 


1.375961 


-0.202161 



Table 1: Monte Carlo and large N data for the lattice step scaling function at uq = 
1.0595. 



For technical details about our algorithm and estimators used in the sim- 
ulations we refer the reader to Ref . [Hj . Results are listed in Table [T] The 
0(3) data for L/a < 64 has been taken from |6., where Seefeld et. al. study 
cutoff effects in E(wq,3, a/L) at the popular value u = 1.0595 already ap- 
pearing in [Hj. We have extended their data to the lattices L/a = 96, 128 
and have added results at iV = 4, 8 to explore the transition to the large N 
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behaviour that we later study semi-analytically. The bottom line of the table 
contains proposed exact continuum values published in [§| for N = 3 and the 
values for N = 4,8 pH]. They depend on some mild theoretical assump- 
tions, but look consistent with our numerical values [HIE]- More precisely, 
an extrapolation with the fit-forms considered below (eqns. (jBJ) - (jHJ) leads 
to extrapolated values that agree with the conjectured exact values within 
errors (propagated from the Monte Carlo statistical errors) as soon as the x 2 
of the fit is acceptable 1 . More details can be found in [Hj. In this publication 
we assume the conjectured continuum values to be correct and use them to 
constrain our extrapolations. We notice that the case N = 8 is very well 
described by the large N expansion. Values of S at N = 8 and N = oo 
differ by about 2.2%. This discrepancy is lowered to about 0.4% by the 1/N 
correction that we shall derive in this paper. 

We have a first look at cutoff effects in Fig. Q where we simply plot E/cr 
versus (a/L) 2 for the different values of N together with the exact curve for 
N = oo. Again N = 8 is already close to the large N limit. The lattice 
artifacts for N > 3 generally turn out to be much smaller than for N = 3, 
which makes the numerical investigation of their structure more difficult. 
While the N = 3 data look curved to the naked eye in this representation, 
this cannot be decided so easily for the higher iV-values. In order to analyze 
the cutoff effects in more detail we will fit several analytic forms to the data. 
In the following Section we explain what kind of fits we have done and how 
and what we can learn from them. 

2.1 Fits to the data 

Symanzik's analysis of the cutoff dependence of lattice Feynman diagrams 
suggests that leading lattice artifacts are quadratic in the lattice spacing 
modified by slowly varying logarithmic functions. It should be noted how- 
ever that logarithms generically appear in an Z'th order polynomial at I loop 
order. If the series were summable the leading power could thus be modi- 
fied. Nevertheless in typical lattice applications it is assumed that the integer 
power behaviour is modified only weakly and that this can be neglected for 
most attainable statistical precisions in an extrapolation. It was first found 
in [H] that in the 0(3) model instead of the expected quadratic dependence a 
behaviour more close to a linear is suggested. With this background we first 

1 There is one exception to this for N = 3: a free exponent fit (form including 
coarser lattices, where the conjectured value is missed by several standard deviations with 
an exponent close to one. 
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Figure 2: Results of the fit to eq. JSJ for N = 3. In the upper half the optimal exponent 
and the limits where the fit with fixed a becomes unlikely is plotted. The lower half shows 
Q(<*opt)- The dotted lines are to guide the eye. 



investigate fits with a general exponent a 

A: c + Cl (a/L) a , (6) 

where Cq is fixed by our exact continuum data. Such a fit results in an optimal 
exponent a opt that minimizes the x 2 deviation between (JBJ) and the data. To 
judge the plausibility of such a fit we compute the probability Q of finding a 
value of x 2 larger or equal to the one encountered in the fit at hand (goodness 
of fit test fTT|). In general, values Q < 0.1 are considered implausible. 

For iV = 3 fit results are shown in Fig.|2]which requires some explanation. 
Since all theoretical forms of the cutoff behavior are only asymptotic, we 
successively exclude coarser lattices from the fit. We therefore find values a 
and Q as functions of (L/a) m i n , the size in lattice units of the coarsest lattice 
still included in the fit. The top labels list the numbers v of the remaining 
degrees of freedom in the fit. Errors of a opt are assessed in two ways. The 
errorbars are limited by those values of a where the optimal x 2 grows by 
unity. The plus signs are the limits where Q falls below 0.1. When there are 
no fits with Q > 0.1 no errors are drawn. We see here that acceptable fits are 
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(L/a) m i n 
Figure 3: As Fig.EJbut for N = 4. 



only achieved with lattices L/a > 24. Moreover they are neither consistent 
with the Symanzik value a = 2 (dash-dotted line) nor with a = 1 but rather 
favor values in between. For N = 4,8 the equivalent plots are shown in Fig. El 
and Fig. HI In these cases the proposed fit is acceptable for all (L/a) mm > 10. 
For N = 4 the preferred a- values are still smaller than two, but N = 8 starts 
to look consistent with the whole range between one and two. 

In addition to fit A with a general correction to scaling exponent we 
also tried the Symanzik and the ad hoc linear form with both including the 
simplest logarithmic correction 

B : c + d (a/L) + c 2 (a/L) \n(a/L) . (7) 
C : c + ci (a/L) 2 + c 2 (a/L) 2 \n(a/L) , (8) 

Note that at iV = oo the form C is known to apply (see next section). In the 
upper and the lower subplot of Fig. El the Q-values are shown for the linear 
(eq. (jZJ) and quadratic (eq. (JBJ)) form respectively. For iV = 4, 8 both fits 
are tolerable in the whole (L/a) m j n -range. In the case N = 3 and at small 
(L/a) m i n both forms become unlikely. For (L/a) min = 16 the form C, with 
quadratic lattice artifacts, gives the better fit, but for (L/a) min = 12 it is 
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(L/a) m i n 
Figure 4: As Fig. but for N 



already intolerable whereas the linear fit with Q m 10 _1 could perhaps still 
be accepted. 

In summary, the single exponent ansatz A hints at 'effective' values be- 
tween one and two in the range of cutoffs studied. With logarithmic mod- 
ifications inluded, linear and quadratic dependence is asymptotically about 
equally well compatible with our data. All fits studied here have at most two 
free parameters. More complicated fits can be extended to coarser lattices, 
and then the additional parameters are actually largely determined by these 
additional included data. We do not think that much can be learnt from 
such more involved fits. We rather turn now to the study of large N to then 
apply the same fitting procedure in a case where we have additional insight 
in the true cutoff dependence. 



3 O(N) model at large N 

At finite N Monte Carlo simulations were our tool to go beyond perturbation 
theory. In this section we complement this information by evaluating the 
first two terms of the 1/N expansion of the step scaling function. As we 
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Figure 5: Q- values of the fits to eq. Q (linear) and eq. JHJl (quadratic) against (L/a) m i n . 
For N — 4,8 both fits are perfect in the whole range. For N = 3 both fits start to become 
unlikely at (L/a) nl - m = 12. 



shall discuss they are nonperturbative in the renormalized coupling M(L)L. 
Related large N computations are found for example in [12] and [THj, but to 
our knowledge the step scaling function as defined here and in particular its 
subleading correction is not in the literature. 

For simplicity we employ lattice units in the following with the lattice 
spacing set to a = 1. Scaling is then studied in 1/L. 



3.1 Spin correlation function 

To extract the mass-gap we need to study the spin two-point function. We 
Fourier represent the spin constraints in (0) and add a source term to obtain 
the generating functional of spin correlation functions 

Z(J) = [ y\ d N s ^L e -S(s)-( a ,s>-l) + (J,s) (9) 

/ - LJ - 2ni 

J X 
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with each a(x) integrated along the imaginary axis. We rescale a, shift 
contours a(x) = ml + iip(x) and perform the Gaussian s-integrals to obtain 



z^oc j n 



dlpe S eS (<p)+^(J,K-ij) (1Q) 

x 

with 



S eff M = ^trlogX-^5> (11) 

X 

involving the operator 

K = -d^d; + ml + i(p = (IT 1 + i<p) = D-\l + iD(f), (12) 

where the scalar field f is regarded as a diagonal matrix and 

D = (-d tl d; + m 2 )- 1 . (13) 

The mass m\ is chosen such that a(x) = m\ is a saddlepoint or equivalently, 
that S e s(ip) is minimal for if — 0. This occurs if the gap equation 

^ = ^tr^U = itrD (14) 

is fulfilled with the volume V equaling the number of lattice sites. 

Next we expand S c s in powers of ip, omit a constant and rescale ip — > 
y/2/N(p to find 



Ses((p) 
Similarly we then get 



\ tr(^) 2 - tr(^) 3 " ^ tr(D^) 4 + ■ ■ ■ ■ 



(15) 



K 1 = D — iy ^D V D - ^D V D V D + .... 
The spin two point function is given by 

(si(x)sj(y)) = jy5ijG{x,y), 

G(x,y) = (K- 1 (x,y)) 

with the expectation value taken with S c g and employing an obvious kernel 
notation for the operator K~ l . Upon 1/N expansion this implies 

G = D - A (DipDipD), + A (ity£> tr(^) 3 > + 0(iV 2 ) (16) 

9 



FIGURE 6: Diagrammatic representation of (|21(l and (|22() . 



with (• • -) taken now with only the quadratic part of S e g. We introduce the 
selfenergy operator H = j^Hi + 0{1/ N 2 ) by setting 

G- 1 = D- l + H (17) 

and read off by comparing with (fTtjj) 

H x = 2 (<pD<p) - (<p tr(Z^) 3 ) = H lta + H 1;b . (18) 
Going to momentum space with 

D(jp) = 1 ^-— 2 , ^ = 2sin(^/2) (19) 
and the (^-propagator VT(p) determined in the first term in |T5|) as 

9 



the two contributions have the form 

2 
V 



= ^1] £>(<?) W%-P) (21) 



and 

#1,6 = -W(0)-^ ^ 1^(9)^(9 - r)D 2 (r). (22) 



V" 2 



These terms are visualized by the diagrams in Fig. |U] where D is represented 
by solid and W by wiggly lines. Using the explicit form of D we can eliminate 
the double sum in the constant Hxjb and find 

K b = W{^Y. W ^^f~ 1 ^)- ( 23 ) 
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3.2 Mass gap 

For our two-dimensional volume V = TL we now take the limit T — > oo. 
The momentum po becomes continuous and we have to replace 

— X^ 1 f dp ° 

V 2^ ■ ■ ■ ■ ~* J 2vr " " 

V Pi 

The gap equation becomes in this limit 

1 V f d ^ 1 - 1 V 1 = 1 (24) 

L A.y 2 7rp 2 + m 2 l L^2^(p 1 ) v /l+^ 2 (p 1 )/4 7 



with o) 2 (pi) = (p 2 + m§). It furnishes a one-to-one relation (at given L) 
between Wq and 7 > 0. 

The mass gap M > of the transfer matrix can be extracted from the 
two point function in momentum space G(p) by finding the smallest real M 
with 

G- 1 (p M ) = 0, p M =(iM,0). (25) 

We expand 

M = M + —Mi + .... (26) 
At leading order G coincides with D which yields 

M = ln l+m 0A /l + m 2 /4 + m^/2 =M {j,L). (27) 



If we now combine (I25j) . (j26j) . (fTTj) we obtain for the first correction to the 
mass 

M x = H ) iPMo) (28) 
2m 0V /l + mg/4 

In the practical evaluation M (7,L) follows easily from solving ()24j) for 
instance by the Newton Raphson algorithm. For the correction it is essential 
that also the T — > 00 limit of W can be taken in closed form. In a first step 
we write it as a contour integral in z = e m over the unit circle 



w-\p) = \Y,j 



dz 



Here we have introduced 



2niz 2 cosh(#) — z — z' 1 2 cosh(<5) — zz 



z 



z = e" ipo , cosh(5) = 1 + w 2 (gi - Pl )/2, cosh(fl) = 1 + Co 2 {qi)/2. 
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After some algebraic simplifications we arrive at 



With this form the computation of both Hi A (p Mo ) and f/^f, in (J2~T|) and (|2*S)l 
leaves us with a single numerical integration over q e [0, 7r] (using symmetry) 
in combination with a two-fold summation over the spatial momenta. This 
is easily feasible on a PC with high precision and up to L of a few hundred. 
Between the two diagrams a divergence oc L 2 cancels 2 . 

A final remark concerns the continuation of H\ a to imaginary momentum. 
For real momenta p we can obviously use the shifted form 

£i,.(p) = t£ / j^D(q + r)W(q-p + r) 

instead of (|21j) with any admissable lattice momentum r (real ro, r\ = 2 r nnjV) 
including the case r = p. Only in the latter case we encounter a singularity 
if then we simply substitute p — > pm and a more careful continuation would 
be required. We used (|21|) in its original form amounting to r = 0. 

3.3 Step scaling function 

We can now determine the coefficients of the expansion (JHJ). For the leading 
term we read off [14| 

V (u,L) = 2LM { 1 ,2L)\ LMo{ ^ L)=u (30) 
with 7 determined by u. For the next term we find (with the same 7) 

E x (ti, L) = 2LM 1 ( 1 , 2L) - LM^j, L) d ^ L) , (31) 



where the derivative is evaluated in the form 

dZ (u,L) 2 



9M (7,2L) 

81 (32) 



Q u " 9Mo( 7 ,L) • 

The derivatives dMo/d'-f follow from ()27)1 and 7-derivatives on both sides of 
()24j) avoiding the need to take any derivative numerically. 

Each term in the expansion is expected to have a continuum limit for 
fixed u 

(Ti(u) = lim T,i(u,L) (33) 

L— >oo 



2 We thank Janos Balog pointing this out. 
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3.4 Infinite N 
3.4.1 Continuum 

At leading order the step scaling function is given just by the gap equation. 
In ref . ^2] we find asymptotic large L expansions of lattice sums relevant for 
that yield 

- = - [InL + G (u/27r)] + - + C- + 0(1/L 2 ). (34) 

7 7T U 07T Z L z 



Here G is some constant and Go is given by 



G Q (a) = 



n=l 



1 



\Jn 2 + a 2 n 



(35) 



Note that the Symanzik form of the leading lattice artifacts holds here for 
all couplings. The exact continuum step scaling function ao(u) is given by 
the transcendental equation 

- [In 2 + Go(<x /27t)] + — = -G (u/2n) + -. (36) 

7T Gq H U 

The function Go possesses a perturbative expansion, convergent 3 for a < 1, 

= E ¥§^ 2n + 1 )(-« 2 /4)" (37) 

n=l ^ ' 

which implies the expansion coefficients of a to all orders. The first ones in 

a = u + sqU 2 + Siu 3 + . . . (38) 

are 

In 2 2 3 4 C(3) 5 2 7C(3) 

s = , s 1 = s , s 2 = s , s 3 = s -s -—, s 4 = s -s — — . (39) 

7T 47T° 07T° 

In ref.^3] the four loop /5-function at finite iV was computed for the LWW 
coupling. Its iV — » oo limit yields an expansion for the step scaling function 
that coincides with the present series and from S3 we find for a constant 
that was determined numerically before the compatible result xi — — C(3) = 
—1.20. . .. Note that for the quantities considered here we confirm the in- 
terchangeability of the large N- and the cutoff-limit leading to renormalized 
perturbation theory. 

3 This is closely related to the series for M(L)/A-^ discussed in |7j. 
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Figure 7: Exact infinite N step scaling function ctq together with perturbative and large u 
approximations. Labels 1,2,... denote the expansion (|38[) including terms up to so, s\, . . . 
and agree with the loop-order in conventional perturbation theory. 



At large a the function (|35|) has an asymptotic expansion 

j oo 

G {a) = - ln(a/2) - — + 2 V K {27ma) + 7e (40) 



2a 

n=l 



with Euler's constant 7e and the modified Bessel function Kq. It implies the 
usual exponential corrections to the finite volume mass gap 



a Q (u 



2u [1 - 2K (u)) =2u\l- ^/2~Kjue^(-u) \ . (41) 



In Fig. Owe plot the exact continuum step scaling function together with its 
small and large u expansions. We remark that a few orders of perturbation 
theory approximate cr well up to u ^ 1.5. 

Eq. (jHUj) can of course also be analyzed for other rescaling factors then 
two. In particular from infinitesimal rescaling we obtain the nonperturbative 
/3-function in the form 



u 2 



f3(u) = ^— - - + 0(1/N). (42) 

71 1 - i?G' {u/27i) 
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3.4.2 Lattice artifacts 

At finite L we compute T, (u,L) exactly by using (J24j) twice (with L and 
2L), eliminating 7 and expressing mo by Mq via the inverse of (|27jl . itlq = 
2 sinh(Mo/2). This is implemented numerically to machine precision. An- 
other possibility is the small volume expansion of the gap equation 
yielding 



2 

7 u 



(43) 



with 



ly \ d 3 = -J-y 1+ ^ /2 (44) 



4 = -6^' 4 =36oW (45) 
Numerically we find by the method of App. D of ref.[TE] 

cfi(L) = — - 0.0703275870013 + 0.0872665^7 + 0(L" 4 ) (46) 

7T L 2 

and 

X In -Z^ 1 
d 3 (L) = -0.0048460224504 - —-yj- + 0.0220539— + 0(L~ 4 ) (47) 

with errors expected to lie beyond the quoted digits. 

In view of applications of the step scaling technique for QCD it is interest- 
ing to compare the exact lattice artifacts with their perturbative expansion. 
We form the relative deviation 

6{u, L) = = 5 x {L)u + 5 2 (L)u 2 + .... (48) 

Setting 

Adi(L) = di{2L) - di(L) (49) 
and combining series we find 

S 1 (L) = Adx(L) - s , (50) 
S 2 (L) = Ad 2 (L) - si + s 2 + Adi(L)5i(L), (51) 
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Figure 8: Deviation of Eg from the continuum limit at uq = 1.0595. 



and similar but lengthy expressions for 5 3j 4(L). In Fig. [HJ the exact S(L) is 
plotted together with perturbative approximations at uq = 1.0595. While 
the continuum limit in this range is well described by perturbation theory, 
this is not really the case for 5. The same was observed in [Hj at iV = 3. 
The leading term does not even have the right sign. This situation improves 
at yet smaller couplings as shown in Fig. El Heuristically, with 5 being a 
short distance quantity, we thought it might be a good idea to reexpand 
it in the bare coupling j(u, L) which is always known nonperturbatively in 
computations of step scaling functions. This experiment did however not 
lead to a more accurate perturbative description of 5 in the present case. 

3.5 1/N correction 

For a number of u- values we numerically determined £i(iz, L) up to L = 128 
and higher for Uq. Larger values are possible with some effort, but are of 
limited use unless also the precision is raised beyond 64 bits. Sofar we have 
not succeeded in analytically extracting the large L-behavior and therefore 
perform a numerical analysis here. 
The expected form is 

^(u,L) = a 1 {u) + ^\u,L) + jjZ { i\u,L) + 0{L- & ) (52) 
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Figure 9: As Fig.® but for u = 0.5. 



where the T\ have an only weak residual L-dependence. For the entirely 
analogous expansion of the leading S it rigorously follows from that 
Sq 2 '* is a linear and £ 4 ^ a quadratic polynomial in ln(L). 
In terms of the symmetric difference 

A/(L) = i(/(L + l)-/(L-l)) (53) 
we define similarly to JTTj 'blocked functions' 

^L 3 AEj (54) 

— L 4 (2 + LA) 2 LAS i (55) 
16 

(2) 

for i — 0, 1. The reasoning behind this construction is that in Bq the contin- 
uum part of So has been canceled and the L~ 2 components are promoted to 
order unity. The normalization is such that the coefficient of the term linear 
in InL in £q 2 ^ is the same as in Bq . Similarly in B^ the L~ 2 components 
are canceled and the coefficient of (lnL) 2 /L 4 is preserved. Although the a 
priori form of is not analytically known sofar, we construct B± in the 
same way In Fig. [TUl we plot Bq at u = u for a range of L and see the 
expected slowly varying L dependence. The analogous plot Fig. [TT] for the 



Bf = - 
Bf = - 
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1/N correction shows a rather similar structure. Both lower plots eventually 
start to become irregular at large L due to roundoff noise. For both terms 
Sp' 1 looks strictly linear in InL and the (InL) 2 components in S^, i. e. the 
curvature at large L, seems to be very weak. For other w-values the situation 
is qualitatively similar, although coefficients are larger for larger u. 

After this more qualitative investigation we obtain the coefficients in the 
asymptotic expansion 

E 1 (u, L) = a x {u) + A(u) 1 ^ + B{u)^- 2 + 0(L" 4 ) (56) 

by the method introduced in [XH|. The resulting numbers are collected in 
TabE] and the continuum values are plotted in Fig. El The correction ob- 
viously has to (and does) vanish for u — > 0, oo since a(u, N) assumes N- 
independent trivial limits u and 2u. 

It seems rather difficult but perhaps not hopeless to analytically derive the 
asymptotic L-dependence of Si. We hope to come back to such an analysis 
in a future publication [TKj . 

It is instructive to perform with our large N data the same analysis as 
for the Monte Carlo data. For this purpose we 'fuzz' our data with artificial 
errors. For instance we add to Ei(wq, L) in the range 10 < L < 128 Gaussian 
random numbers of width 2 x 10 -5 . This size is chosen such that the lattice 
artefact |Si(m , 10) — ai(uo)\ of the 1/N correction gets a relative 'statistical' 
error similar to the N = 3 case. The resulting plot analogous to Fig. |21 is 
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FIGURE 11: As Fig. E3 but for the subleading term Si 

shown in Fig. We see a rather similar picture for this case where we 
are confident of the lattice artefacts really having the Symanzik form. Very 
similar plots can also be generated from the leading N = oo term. Also in 
Fig. |U the new data points would resemble the structure exhibited by the 
Monte Carlo data. 

4 Conclusions 

In this work we have further investigated the scaling violations of the finite 
volume mass gap step scaling function of the O(N) model at various N. 
The discretizations are restricted to the standard action. The behaviour 
at N = 3 is confirmed by our added data and by our method of analysis. 
With fits of the form © a rather large number of coarser lattices has to be 
excluded, and then a cutoff dependence with an exponent 1 < a < 2 seems to 
be suggested in the accessible range of lattice spacings. With one additional 
logarithmic correction both exponents 1 (ad hoc) and 2 (Symanzik) (|7|8jl can 
be used with again only the finer lattices included. For N = 4, 8 a similar but 
less pronounced picture emerges. Our sensitivity is weaker since the lattice 
artefacts are generally much smaller and hence harder to measure. Although 
there may be a trend of 'normalization' with larger N, also at N = 4 our 
data point toward a < 2 unless all lattices with L/a < 24 are excluded. 
It would be very interesting if theoretical principles could be found to 
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u 


<y\ 


.4 


B 


0.5 


-0.055886 


0.012 


-0.066 


1 


-0.187974 


0.220 


-0.184 


1.0595 


-0.202161 


0.284 


-0.200 


1.5 


-0.227820 


1.27 


-0.442 


1.75 


-0.158761 


2.23 


-0.84 


2 


-0.046463 


3.26 


-1.56 


2.25 


0.072185 


4.13 


-2.52 


2.5 


0.167824 


4.67 


-3.53(1) 


3 


0.258377 


4.82 


-5.12(1) 


3.5 


0.251671(1) 


4.21(1) 


-5.71(3) 


4 


0.205320 


3.38(1) 


-5.50(4) 


4.5 


0.152853 


2.58(1) 


-4.85(4) 


5 


0.107791 


1.92(1) 


-4.04(4) 


5.5 


0.073422 


1.40(1) 


-3.24(4) 


6 


0.048858 


1.01(1) 


-2.52(4) 



Table 2: Continuum extrapolation of Si, see eg. 1(55)1 



organize a leading log summation of leading cutoff effects in analogy to the 
renormalization group improvement of continuum results. Such an attempt 
is here left to future study. 

At iV = oo it is rigorously known due to [12] that cutoff effects are of 
the Symanzik type with just 1/L 2 , ln(L)/L 2 contributions. The leading 1/N 
correction to the step scaling function, worked out here for the first time, 
somewhat to our surprise seems to have just the same behaviour, although 
this can be only demonstrated numerically at present. Needless to say, we 
cannot exclude other admixtures if they are only small enough. When these 
exact results are analyzed in the same way as the Monte Carlo data, a qualita- 
tively similar picture emerges. It is hence very clear that we cannot conclude 
a truly anomalous scaling behavior from the latter. 

The large iV sigma model proved to be a very interesting case, also because 
we here have an analytic handle on quantities that are nonperturbative in the 
coupling of an asymptotically free theory. This is true for both the continuum 
and at finite lattice spacing. This simplified laboratory can certainly be useful 
to study other issues. 
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FIGURE 12: Leading correction <7i(u). 
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FIGURE 13: Fit A for the large N data Si(uo, L) with artificial errors added. 
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